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ABSTRACT 

We consider the evolution of the sky-averaged 21 cm background during the early phases 
of structure formation. Using simple analytic models, we calculate the thermal and ionization 
histories, assuming that stellar photons dominate the radiation background. The resulting 21 
cm spectra can constrain the properties of the first generations of stars and quasars. If Popula- 
tion II stars dominate, Lya coupling renders the IGM visible before it is heated by X-rays and 
long before reionization. Thus the 21 cm background has a strong absorption epoch followed 
by weaker emission that fades during reionization. The harder spectra of very massive Popu- 
lation III stars compress these transitions into a shorter time interval and decreases the signal 
amplitude. However, the reionization epoch will remain visible except in extreme cases. The 
global 2 1 cm signal will be challenging to observe because of astronomical foregrounds, but it 
offers an exciting opportunity to study the first sources of light. It also fixes the overall ampli- 
tude of the fluctuating background whose detection is a major goal of several next-generation 
low-frequency radio interferometers. 
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1 INTRODUCTION 

■ The earliest generations of cosmic structures have profound effects 
on the Universe around them and on later generations of stars and 

, galaxies. These sources "reionize" the intergalactic medium (IGM), 

' enrich galaxies (and the IGM) with the first heavy elements, and 
seed larger galaxies. 

Observations now paint a fascinating - but confusing - 
picture of their properties. Unfortunately, observing "typical" 

, galaxies at z > 8 is beyond th e capabilities of existing instruments 
(though see iBouwens et alJ |2004), so current constraints use 
reionization as an indirect proxy. The rapid evolution of the mean 
Lya optical depth in quasar absorption spectra seems to indicate 
that reionization is ending at 2 ~ 6 JFanet alJl200lL 120051) . 
Corroborating evidence comes from the proximity zones around 
these quasars, whic h also imply rapid evolution in the mean 
ionized fraction Xi 1 Wvithe & Loeb '20041 iMesinger & HaimanI 
[2004; Wvithe etal. 2005; Fan et al. 200^ as well as large-scale 
fluctuations in the transmission. But both are difficult to interpret: 
the latter, for example , may also be ch aracteristic of reionization 
iWvithe & Loebll2005t [pan et alj|2005l) . but so far the variation s 
are also consistent with density fluctuations JLidz et all l200d) . 
On the other hand, large-scale polarization of the cosmic mi- 
crowave background (CMB) observed by WMAP implies an early 
beginning to reionization at z > 10, al beit with large uncertainties 
(Page et al. 2006; Sperg el et alj|200a) . How these puzzle pieces 
fit together - and, more fundamentally, how luminous sources 
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evolve, how feedback processes (such as metal enrichment, 
photoionization, or other radiation backgrounds) influence them, 
and h ow small-scale structure evolves in the IGM- remain myste- 
rious <Wvithe & Loebll2003l ICerJl2003lJ; IWvithe & Loebll2003bl: 



"Cen '20033' 



Tl 



Haiman & Holder 2003; Somerville & Livio 2003 
Pukugita & Kawasaki 2003 ; Onken & Miralda-Pscude 2004 



Purla netto & Loe b 2005; Choudhurv & Perraral 1200^ llhev et all 
■2005D . Unfortunately, identifying the relevant feedback processes 
and uncovering the source properties are impossible with existing 
data. 

One promising method to disentangle these degeneracies is 
the 21 cm transition of neutral hydrogen (see lPurlanetto et all200a 
for a review). Before reionization, the neutral IGM gas should 
be visible against the CMB through this transit ion iPield 195 3 
ISunvaev & Zeldov ich"l972";"Hogan & Reesll979l) . In principle, we 
can use i t to map out the growth of the cosmic web and HII regions 
at z > 6 JScott & ReeJll990tlMadau et alJll997l;IZaldarriaga eTall 
I2OOJ). The key advantage of the 21 cm transition is that it is a line: 
hence every observed frequency corresponds to a different distance, 
and we can isolate each phase in the IGM's history. 

Of the many applications of the 21 cm line, the simplest (at 
least conceptually) is to measure the global background as a func- 
tion of frequency. This would allow us to measure Xi (z) and, at ear- 
lier redshifts, the growth of soft-ultraviolet (UV) and X-ray radia- 
tion backgrounds, which are difficult to constrain in any other way. 
Three quantities are relevant: the thermal history, the Lya back- 
ground (both of which determine the spin temperature Ts of the 
IGM), and the ionization history (which determines when the signal 
fades). Unfortunately, this global background is likely to be diffi- 
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cult to measure because of contamination by Galactic synchrotron 
radiation (see ^. As a result, it has received little theoretical at- 
tention, and most existing work has simply argued that Ts 2> 
T-y (where T~, is the CMB temperature) during reionization (e.g., 
i Madau et alJl997l:ICiardi & Madaul200a:IChen & Miralda-Escudd 
|2004) . This is unfortunate because a number of efforts are never- 
theless underway to observe it (and the challenges are not neces- 
sarily much greater than observing fluctuations in the background). 
Indirect techniques may also allow us to measure the background 
JBarkana & Loebl2005bl:ICooravl26o3) . And, of course, the global 
signal sets the overall amplitude of fluctuations, so those should 
be considered within this context. One exception is Sethi ( 2005), 
who studied the evolution of the global background in a limited set 
of models that did not fully explore the many uncertainties in the 
properties of the first sources. 

The purpose of this work is to compute the high-redshift 21 
cm signal and to identify those properties that can be measured 
from the mean background. Because we know so little about the 
first sources, it is of course impossible to make robust predictions 
for the background. We will instead parameterize the relevant pro- 
cesses and consider how they affect the mean signal. We first give 
some background on the 21 cm transition in fJ2| We then describe 
processes relevant to the thermal history of the IGM in fj3| to the 
spin temperature in fJ4| and to the ionization history in ij5| We then 
describe some general consequences for the 21 cm history in S|6| 
and examine some specific models in [Q We conclude in SJ8| 

In our numerical calculations, we assume a cosmology with 
n™ = 0.26, Q.K = 0.74, Q,b = 0.044, H = lOO/i km s"^ Mpc"^ 
(with h — 0.74), n = 0. 95, and erg = 0.8, consistent with the 
most recent measurements iSpergel et al. 2006), although we have 
increased as from the best-fit WMAP value in order to improve 
agreement with weak-lensing data. 



s — 1 ~ ■ 

_L -p Xa T" Xc 



(2) 



Here Xe oc n is the colli sional coupling coeffic ient for H-H inter- 
action s iZvgelmaril2005l) and H-e~ collisions iSmithll 1 966t iLisztl 
[2001). The middle term describes the Wouthuysen-Field effect, in 
which absorption and re-emission of Ly o; photons mixes the hy- 
perfine states. The coupling coefficient is iChen & Miralda-Escuda 
l2004lHiratj2005h 

Xc ^ 1.81 X 10^^(1 + z)-^ScJc, (3) 

where Sa is a factor of order unity describing the detailed 
atomic physics of the scattering process and Ja is the back- 
ground flux at the Lya frequency in units cm~^ s~^ Hz"^ sr~^; 
the Wouthuysen-Field effect becomes efficient when there is about 
one photon per baryon near this frequency. It couples Ts to an ef- 
fective color temperature Tc (in most circumstances, Tc ~ Tk', 
^^ield 1959). We use the numerical fits of Hirata 1 20051) for Sa and 



3 THE THERMAL HISTORY 

To compute the evolution of 5Tb, we therefore need Tk (z), Ja {z), 
and Xi{z). In the next three sections, we will examine each in turn, 
beginning with the kinetic temperature. The evolution equation is 

'J^ = -2Hiz)TK+^-J2^, 
at S ^ — ' kbu 



(4) 



where the first term on the right hand side is the p dV work from 
the expanding universe and e^ is the energy input into the gas per 
second per unit (physical) volume through process i. We include a 
number of possible input mechanisms. 



2 THE 21 CM TRANSITION 

We review the relevant characteristi cs of the 21 cm tran sition here; 
we refer the interested reader to Fu rlanetto et alj yOOg) for a more 
comprehensive discussion. The 21 cm brightness temperature (rel- 
ative to the CMB) of a patch of the IGM is 



ST, 



27xm{l + 5) 
V Ts 



0.023 



)( 



0.15 l + z\'^^^ 



n^K^ 10 



mK, 



(1) 



where 5 is the fractional overdensity, xhi = 1 — x^ is the neu- 
tral fraction, and Xi is the ionized fraction. Note that the patch 
will appear in absorption if Ts < Tj and emission otherwise. 
Here we have assumed that the patch expands uniformly with the 
Hubble flow; radial peculiar velocities also affect ST, by changing 
the mapping from distance to frequency ( Bharadwai & Ali 2004 ; 
iBarkana & Loebl 12005 ah . However, these anisotropic fluctuations 
do not affect the mean signal, so we will ignore them. We will use 
Xi to denote the globally averaged ionized fraction; the tempera- 
tures are likely to be much more uniform at most times, and we 
will implicitly treat the m as independent of position (though see 
iBarkana & Loebl2005hl) . 

The spin temperature Ts is determined by competition be- 
tween three processes: sc attering of CMB photons, collis ions, and 
scattering of Lya photons iWouthuvsenI 1 95 2l ; lFielci 1 95 8h . In equi- 
librium iHirata 2005i) , 



3.1 Compton Heating 

At high redshifts, Compton scattering between CMB photons and 
the resi dual free e lectrons in the IGM dominates. The heating rate 
is re.g.. lSeageret al. 2000) 



2 e. 



Xi 8aTU-y 

3 ksn 1 + /ho + Xi SrrieC 



{T-,-Tk), 



(5) 



where /hc is the helium fraction (by number), u~j oc T^ is the en- 
ergy density of the CMB, and or is the Thomson cross-section. The 
first factor occurs because CMB photons scatter off free electrons, 
but the heat must be shared with all the particles. The efficiency of 
Compton heating decreases with cosmic time as electrons recom- 
bine and u^, falls; the IGM thermally decouples from the CMB at 
(Peebles 1993) 

1 + 2dec ~ 150(!:^b/^V0•023)^/^ (6) 

Below this point, Tk oc (1 + z)^, as appropriate for a thermally 
isolated expanding gas. 



3.2 X-ray Heating 

The thermal evolution becomes much more complex (and uncer- 
tain) once the first luminous sources appear. X-rays from galax- 
ies and quasars likely provide most of the heat for the low-density 
IGM. Given our large uncertainties about the nature of high- 
redshift objects, we will take a simple parameterized approach cal- 
ibrated to nearby starburst galaxies. These objects show a clear 
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correlation between star formation rate (as measured by the to- 
tal infrared luminosity) and hard X-ray (2-10 keV) luminosity 
(Grimm et al. 2003 ; Ranalli et al. 2003 ; Gilfanov et al. 2004D : 



Lx = 3.4 X 10*° /js 



SFR 



IM 



©yr 



ergs 



(7) 



where fx is a correction factor accounting for the (unknown) dif- 
ferences at high redshifts. Note that we have transformed to the 
total (> 0.2 keV) X-ray luminosity assu ming a spectral inde x 
a = 1.5 (where L^ ex u^") in this range I Repha eli et aljri993) . 
This is only an average index, of course, but errors can be subsumed 
into the unknown normalization fx ■ 

There is of course no rigorous justification for extrapolating 
this relationship to high redshifts, but it is certainly reasonable. 
Massive X-ray binaries produce the bulk of high-energy emis- 
sion in nearby starbursts. Such systems form when the first mas- 
sive stars die, only a few million years after star formation com- 
mences. So they should also be present in high-redshift galaxies 
(Glover & Brand 2003). Of course, their abundance must change 
with the metallicity and initial mass function (IMF): if, for exam- 
ple, the IMF is top-heavy, fx could be significantly larger than 
unity. But for the moment we can only speculate about the magni- 
tude of such evolutionary effects. 

A second source is inverse-Compton scattering from relativis- 
tic electrons accelerated in the supernova remnants. In the nearby 
universe, only powerful starbursts have strong enough radiation 
fields for this to be competitive with synchrotron cooling; however, 
the CMB has u^. oc (1 + z)'*, so inverse-Compton emission may 
dominate over synchrotron emission in all star-forming environ- 
ments at high redshift I Oh 2001). We can make a simple estimate 
of the extra luminosity contributed by these fast electrons. Suppose 
that a fraction te of the total supernova energy _Esn is used to ac- 
celerate electrons, and that a fraction /e ~ 0.5 of this energy is 
emitted in X-ray photons (this fraction is determined by cooling 
within the shock; iFurlanetto & Loebl2004) . We also let fsN be the 
number of supemovae per unit mass of star formation. Then 



rSN 



1.6 X 10*" /, 



SFR 



i^SN 



-EsN 

005 0.01 Mq-1 1051 erg 



1 M© yr"^ 



erg s 



(8) 



where we have substituted e^ ~ 0.05 as appears to be appropriate 
for nearby supernova remnants (e.g., Kovama et al. 1995). Thus, if 
a large fraction of the electron energy is emitted in X-ray photons, 
inverse-Compton emission from fast electrons can contribute a non- 
negligible - though not necessarily dominant - fraction of the total 
emission. Conveniently, fx ~ 1 still seems reasonable. 

Extremely massive ( >300 M©) Pop III stars^ could have a 
wildly different normalization, of course: each pair-instability su- 
pernova r eleases ~ 10^^-10 ^*^ erg, much more than normal su- 
pemovae JHeger & WoosleM 12002*). The fraction of X-ray bina- 
ries may also be significantly larger, or the remnants may form 
"miniquasars." Moreover, these stars also produce soft X-rays di- 
rectly, because they have surface temperatures Tcff ~ 10^ K 
JBromm et alJl200ll ; ISchaeredl2002h . To a first approximation, the 



stellar spectra are simple blackbodies, in which case ~ 0.3% of the 
total stellar luminosity is emitted in photons with _E > 0.1 keV. 
This would imply fx 3> 1. However, photons near ~ 0.1 keV are 
absorbed relatively qui ckly, and they affect only the region within 
< 0.5 Mpc of galaxies JSethl2005l;ICen.200 6 ) . Thus it is the harder 
photons with which we are primarily concerned, and these come 
from nonthermal sources: at Tc« = 10^ K, only ~ 2 x 10~^ of the 
Pop III stellar luminosity lies in photons with E > 0.2 keV. 

Because X-ray heating is thus intrinsically non-uniform, we 
must develop some approximation for our global treatment - in 
the simplest terms, we must specify the photon energy range that 
we assume can penetrate into the diffuse IGM. Equation Q in- 
cluded all photons with E > 0.2 keV. In actuality, all photons with 
E < Sthick = 1.5^^(^f( l + 2)/10]^/^ keV are absorbed within a 
Hubble length (Oh 2001), so photons near the bottom edge of our 
range will be absorbed over a relatively small - but by no means 
negligible - volume. By including soft photons, we are attempting 
to balance the highly-penetrating hard X-rays with the more local- 
ized soft X-rays that likely carry more total ener gy. Of course, this 
fluctuating component (e.g., Madau et al. 2004) will lead to non- 
uniform heating ( Sethi 200 5). We defer discussion of these fluctu- 
ations to future work; our concern here is with understanding the 
zeroth-order global evolution. Fortunately, uncertainty in the over- 
all effect of inhomogeneity can be subsumed into our much larger 
uncertainty in fx so long as the spectrum is a reasonably shallow 
power law (as appropriate for nonthermal sources). For example, 
iGlover & B rand ( 20031) used a method much like our own but took 
fx ~ 0.2 because they calibrated to the luminosity in the 2-10 
keV band, which implicitly assumes that soft photons are absorbed 
inside or near their source galaxies. 

X-rays deposit energy in the IGM by photoionizing hydro- 
gen and helium; the hot "primary" electron then distributes its en- 
ergy through three channels: (1) collisional ionization, producing 
more secondary electrons, (2) collisional excitation of He, which 
produces a photon capable of ionizing H, and (3) Coulomb colli- 
sions with thermal electrons. The relative cross-sections of these 
processes determine the fraction of X-ray energy going to heating 
(fx,h) and ionization (/x,ion); clearly they depend on both Xi and 
the initial photon energy. For o ur calculations, we will use th e high 
energy limit (_E > 0.1 keV) of lShuU & van Steenberjjl985h . who 
fit their results to 



fx,h 
fx,ion 



Cl[l-(1 
C2(l 



- X 



y.°-2\b2 



ai\i)ii 



(9) 



Here Ci = 0.9971, ai = 0.2663, &i = 1.3163, C2 = 0.3908, 
a2 = 0.4092, and 62 = 1.7592; the fits are accurate to <2% 
in the range of interest. Although not as good for softer photons, 
the qualitative behavior is similar: when Xi is small, most of the 
energy is deposited in ionizations (fx,h ~ 0.2 for E — 3 keV and 
Xi = 10^'^), but fx,h ^ 1 as Xi ^ 1, because collisions with free 
electrons become more important.^ 

We are now in a position to estimate the heating rate from 
X-rays in the early universe. We first assume that the star for- 
mation rate is proportional to the rate at which matter collapses 
into galaxies, which we assume correspond to dark matter halos 



1 Note that we will use "Pop III" to refer exclusively to this very massive 
case. If Pop III stars have more or less normal initial mass functions, they 
produce histories close to our Pop II case. We use the massive stars only for 
contrast with the "standard" histories and do not mean to pass judgment on 
the relative merits of the different initial mass functions. 



^ Note that here Xi refers to the ionized fraction in the predominantly 
neutral gas (outside of HII regions) whose temperature we care about. In 
the models below, we track this separately from the usual ionized fraction, 
which instead describes the filling factor of HII regions, by assuming that 
only X-rays ionize this mostly neutral gas. 
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with m > rrimin, where rrimi-n is a (redshift-dependent) mini- 
mum mass. We will typically take m-mi-n ~ rn4, where 7714 has a 
virial temperature Tyir = 10 K (the threshold for atomic cooling; 
ISarkana & LoeblbOOlh . We define /coii to be t he fraction of mat- 
ter in collapsed halos with m > nimin', for the lPress & Schechtea 
il974l mass function, it takes the simple analytic form 



/coll = erfc 



Sc{z) 



V2cr(mmin) 



(10) 



where Sc is the linear overdensity at virialization and a^ (m) is the 
variance of the density field when smoothed on scale m. By equa- 
tion 0, the X-ray emissivity is proportional to the star formation 
rate and hence to dfcou/dt; the result is 



ejf 



3 kBnH{z) 



lO^K/x (AZi^%Z^i_^), (U) 



0.1 0.2 0.01 



10 



where /* is the star formation efficiency. We immediately see that 
X-ray heating should be quite rapid, provided that /x ^ 0.1. 

Of course, our assumption that ex oc dfcou/dt may be wrong 
because (for example) some exotic process produces X-rays, or be- 
cause quasars are common, or simply because /* is not a constant. 
However, unless the process of interest is uncorrelated with galaxy 
formation, it will not differ dramatically from that we have dis- 
cussed. On the other hand, quasars could strongly increase the nor- 
malization of fx ■ We have chosen to ignore this possibility because 
extrapolating the known population of bright quasars to higher red- 
shifts yields an X-ray backgr ound with negligible effects on the 
IGM iVenkatesan et alj|200ir) . A number of models posit an en- 
tirely new population of "miniquasars" (e.g., iMadau et alj|2004t 
[Ricotti et al. 2005), but we will conservatively ignore them. We 
briefly discuss how such miniquasars could affect our conclusions 
inJl 



3.3 Lyof Heating 

As photons redshift into the Lyman-series resonances, they scatter 
off of neutral hydrogen atoms in the IGM and deposit energy in 
the gas through atomic recoil. Wh ile ini tially thought to be an im- 
portant he at source ( Madau et alJll997» . Chen & Miralda-Escude 
i2004 ) and'Rvbicki ( 2006) showed that the energy deposition rate 
in repeated Lya scatterings is actually extremely small because 
the photons qui ckly e stablish equilibrium with the gas (see also 
iFurlanetto et alj|2006l for a discussion of the approach to equi- 
librium). The heating rate from higher Lyman resonances is also 
small because those photons are quickly transformed into either 
Lyg photons or continuum p hotons that easily escape to infinity 
iPritchard & Furlanettd 120051) . We therefore neglect this mecha- 



3.4 Shock Heating 

Finally, we must consider heating from shocks in the IGM. At low 
redshifts, shocks play an integral role i n its temperature structure 
jCen & Ostriker 1999; Dave et all200ll) . The shock networks sur- 
rounding filaments and halos can be quite complex ( Miniati et al. 
[2000 ; Kang et al. 2005), but their characteristic properties are easy 
to estimate. The typical shock that has just collapsed at redshift z 
will have a peculiar velocity equal to the total distance it has col- 
lapsed divided by the age of the Universe, «sh ~ H(z)Ri^i(z), 
where Rni is the (physical) nonlinear length scale ICen & Ostriken 



19991) . For a strong shock in a monatomic gas, the corresponding 
postshock temperature is 



T, 



16fes 



16fci3 



H'{z)RUz) 



(12) 



where fi is the mean molecular we ight in atomic units. T his yields 
Tsh ~ 10^ K at the present day JCen & Ostrikeilll99 9^. compa- 
rable to the results of numerical simulations. On the other hand 
Tsh ~ 10* K at z ^ 3, close to the mean temperature of a 
photoionized medium JFurlanetto & Loebll2004h . Thus large-scale 
structure shocks decrease in importance at moderate redshifts. 

However, before reionization they are also important, be- 
cause without photoheating Tk can be small and the IGM sound 
speed may be much less than v^h- Unfortunately, studying this 
regime requires extremely high resolution numerical simulations, 
and the current generatio n appears to disagree on their significance. 
'Gnedin & Shaver (2004) found IGM shocks to be the most impor- 
tant factor in determining STf, at 2: < 15, but other simulations sug- 
gest that filament shocks have a relatively modest effect (compara- 
ble to the signal from minihalos; lKuhlen et all200a ; IShapiro et alJ 
2005). One possible reason for the discrepancy is that the latter two 
studies assumed that only overdense gas could be observed against 
the CMB and so only in cluded shocks bounding the cosmic web. 
iGnedin & Shaved i2004) . on the other hand, assumed that the en- 
tire IGM was visible because of Lya coupling. Thus it is possible 
that the tiny sound speed in the IGM allows a population of weak 
shocks to exist outside the cosmic web. Unfortunately, weak shocks 
in the low density IGM may be difficult to resolve properly in simu- 
lations. Another alternative is that the differences occurred because 
the studies focused on different redshift regimes. 

For now, we will assume the latter and use the simple model 
of IFurlanetto & Loebl i2004) to estimate the effect of shocks i n 
the cosmic web (see alsoJ Nath & Silkl200lLlValageas et al.l2002h . 
This prescription assumes that shocks form around a region once 
its mean linearized density exceeds 5sh = 1.06. This is equiv- 
alent to assuming that they appear when spherical perturbations 
"turnaround" and break off from the Hubble flow (i.e., when the 
flow begins to converge). This spherical picture is approximate, but 
it provides a simple criterion to estimate the fraction of mass that 
has been shocked on its way to halo formation. Conveniently, at 
turnaround, the physical velocity is of course zero, so the peculiar 
velocity is almost precisely Wgh ~ H{z)R^i{z). Thus it reproduces 
the characteristic temperature of z = shocks. In the same way as 
the Press & Schechter ( 1974|) mass function, it also yields the dis- 
tribution of shock temperatures at any epoch. 

Thus, by analogy to equation <10> . the fraction of gas that has 
been shocked to a temperature T > Tsh is 



/sh(> Tsh) = erfc 



5s\i(z) 



^o-(r,h) 



(13) 



where Tsh is related to the mass of the underlying perturbation via 
equation <12> . For example, the mass corresponding to a shock 
with Tsh = T-i is M-y ^ 4.3 x 10^ Mq, which is some- 
what greater than the minimum collapse mass (the so-called "filter 
mass;" Gnedin & Hui 1998) - though note that Ssh < Sc so more 
gas will be in the shocked phase than in collapsed objects of a com- 
parable mass. For our model, we will simply compute the fraction 
/sh of gas that has been heated to Tk > T^. We will assume that 
this gas emits strongly relative to the CMB (Tg ^ T-y), while the 
rest remains at the mean temperature (without shocks). Note that, 
once X-rays have heated the IGM above T^, shock heating will 
become less significant than our prescription assumes; however, in 
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that regime the shocked gas and background IGM emit at nearly 
the same level and so shocks have almost no observable effect. 

We emphasize that this treatment is meant to include only 
shocks that are part of structure formation (sheets, filaments, and 
halos). As described above, there may be a population of weaker 
shocks in the low-density IGM that deposit a significant amount of 
thermal energy there. In the language of our model, these would 
correspond to shocks that appear before turnaround due to random 
IGM flows. Our model may thus underestimate Tk at early times 
and overestimate the amplitude and duration of absorption. But, as 
we will see. X-ray heating generally becomes significant at z > 15, 
when such perturbations would presumably be significantly smaller 
than those of iGnedin & Shaver: t2004V 



4 THE SPIN TEMPERATURE 

Once Tk is known, the next step is to co mpute the spin te mpera- 
ture of the gas. We include H-H collisions ( Zyeelman 2005), H-e^ 
collisionsjSmithl26^^Lisz^0^^, and the Wout huysen-Field ef- 
fect iWouthuvseiJI 1 952t iFieldl 19^ : lHiratall2005h as described in 
SJ2| Collisions are unimportant for the low-density IGM during the 
epochs of interest,"^, so the Wouthuysen-Field mechanism is there- 
fore the most important factor. Its strength depends on the soft-UV 
radiation background, for which we will construct a simple, easily 
parameterized model. We again assume that young stars dominate 
the flux, so that we can write the comoving emissivity at frequency 
V is 



e{u, 



(■ -0 / \ rf/coll 



(14) 



The function etiv) is the number of photons produced in the fre- 
quency interval v ± du/2 per baryon incorporated into stars. For 
spin temperature coupling, we are interested only in photons be- 
tween Lyg and the Lym an-limit.^ We use the spectral models of 
ISarkana & Loebl |2005g), which fit a power law between each 
Lyman-series resonance. For Population III stars, a total of A'^^ = 
4800 photons are produced per baryon in the interval (i^a, /^ll), 
with a relatively flat spectral distribution I Bro mm et all200ll) . For 
continuously forming low-metallicity Population II stars. No, = 
9690 photons are produced pe r baryon, with the spe ctrum falling 
steeply past the Ly/? resonance JLeitherer et alll99Sl) . 
The average Lya background is then 



n=2 
= 2_^ ./rcc(n) / dz 



■w4'^n,Z ),(15) 



47r H{z 



•^ Unless, tha t is, X-ray heating is extremely strong and renders Xi > 0.03 
lNusseJ2005l) . '" 

^ Here we ignore Lyman-Iine photons produced inside galaxies as well as 
line photons produced following IGM ionizations. The former is a reason- 
able assumptions because these photons rapidly redshift out of resonance 
and hence only affect a small fraction of the IGM (much of which is pre- 
dominantly ionized). The second is an excellent approximation for Pop II 
stars but less good for Pop III stai's, because of their hai'd spectra. Never- 
theless, such photons also rapidly redshift out of resonance near the edge of 
the HII regions. 



where Ja is the background from photons that originally redshift 
into the Lyn resonance, !/'„ is the frequency at redshift z' that red- 
shifts into that resonance at redshift z, 



(1 + ^0 



(16) 



{n + iy 



(17) 



and 

1 + 4 _ 

l + z 

is the largest redshift from which a photon can redshift into it. The 
factor /rec(H) Is the fraction of the resulting Lyn photons that actu- 
ally cascade through Lyg an d induce strong coupling ( Hirata 2005 
IPritcha rd & Furlanetto 2005). We truncate the sum at rimax = 23 
(following Barkana & Loeb 2005b) . although the closely-spaced 
n 3> 1 levels contribute only a small fraction of the total flux. 
Finally, we compute the coupling coefficient Xa from equation J3}- 
For estimation purposes, we can approximate equation <15> as 



Ja ~ —- free /* rit A/coll "T— (1 

47r Au 



(18) 



Here /roc is the average probability that a photon in the interval 
{uct, J^ll) is converted into a Lyg photon. It is free = 0.63 and 
0.72 for our Pop III and Pop II spectra JPritchard & Furlanettd 
2005) . We have approximated the spectrum as constant (so that 
£6 ~ Na/Au, where Au = i^ll — J^a) and let A/coii ^ /coii be 
the spectrum-weighted collapse fraction over the appropriate red- 
shift intervals in equation <15> . 



5 THE IONIZATION HISTORY 

The next step in understanding the global history of the 21 cm sig- 
nal is Xi{z). This has been a subject of extensive study over the 
past three d ecades (see, e.g., Barkana & Loeb 2001; Haiman 2004l 
ICiardi & Fe rrara 2005 for recent reviews). As such, we will not at- 
tempt to explore all of its aspects. We will again examine a simple 
model encapsulating the major features in the global evolution. 

We associate reionization with the star formation rate in a sim- 
ilar way to the Lyg and X-ray backgrounds (see eqs. 1111 and 1141 ). 
We write 



~dt " ^''^'*~if^ ^ aAC{z,Xi)xi{z)ne{z), 



(19) 



where (" is an ionizing efficiency parameter, a a ~ 4.2 x 
10~ cm s^ is the recombination coefficient at T = 10 K,'' 
C = (n^) / (rie) is the clumping factor, fie is the average elec- 
tron density in ionized regions, and ^ is an ionizing efficiency pa- 
rameter. The first term on the right hand side is the source term for 
ionizing photons and the second term is a sink describing recombi- 
nations. 

We must now specify (^, mmin (for feou), and C. The ionizing 
efficiency can be written 



C = Anef*feseNio 



(20) 



where /esc is the fraction of ionizing photons that escape their host 
galaxy into the IGM, A^ion is the number of ionizing photons per 
baryon produced in stars, and Aue is a correction factor for he- 
lium. With the exception of ^hc, all of these parameters are highly 



^ Note that we use the case-A value, which essentially amounts to as- 
suming that ionizing photons are absorbed inside dense, neutral systems 
iMiralda-Escude et al. ,200a) . 
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uncertain. The star formation efficiency depends on the dynamics 
of star formation in young galaxies, /esc depends on the distribu- 
tion of gas, dust, and stars within galaxies, and Mon depends on 
the stellar IMF and metallicity. We will therefore content ourselves 
with representative v alues of these paramet ers in our calculations 
(see the discussion in lFurlanetto et all2006l for more information). 

Of course, both C, and mmin could evolve with time be- 
cause of feedback. For example, metal enrichment could decrease 
C, if Pop II stars are less efficient than their Pop III progenitors, 
or nimin could increase becaus e of photoheating during reion- 
izationjkees 1986; Efst athioJ [l992 : Thoul & Weinberg 1 1 996t 
iKitavama & Ik euchi 200C; but see Diikstraet al. 2004). Feedback 
is one promising method to slow reionization and reconcile the 
SDSS quasar data with WMAP, but the mechanisms at play and 
their interactions are essentially unconstrained. For our purposes, 
the important question is how they affect the global reionization 
history. Early efforts suggested that they might intro duce sharp fea- 
tures s uch as "double reionization" (■Cea.2003b. : .Wvithe & Loebl 
12003 ah . However, this requires a sharp drop in the ionizing effi- 
ciency that is synchronized across the entire Universe. That is dif- 
ficult to arrange when the inhomogeneities inherent to local feed- 
back mechanisms are included ( Fur lanetto & Loetl2005l) . Instead, 
feedback is likely to prolong reionization or introduce a "stalling" 
phase. W e will show some example s generated using the feedback 
models of iFurlanetto & Loebl y003) to illustrate its effects. 

The clumping factor can, in principle, be computed through 
numerical simulations with radiative transfer. Unfortunately, that 
requires fully resolving small-scale structure in the IGM as well 
as its evolution throughout reionization. Moreover, C depends on 
how the IGM is ionized - for instance, whether low-density gas is 
ionized first or whether man y photons are consumed ionizing dense 
blobs when Xi is still small jFurlanetto & Ohl2005l) . Thus C must 
be recalculated for each different set of model parameters. Unfor- 
tunately, simulations have not yet achieved these goals. Following 
our usual policy, we will therefore use a simple analytic model 
JMiralda-Escude et al._200 0): in this prescription, voids are ionized 
first. This probably underestimates the mean recombination rate, so 
reionization may be a bit slower than our models predict (and begin 
a bit later). On the other hand, simulations suggest that recombina- 
tions do not dramatically slow reionization, even if minihalos are 
included (Ciardietal. 2005). Note that when recombinations are 
slow, equation <19> can be approximated as 



C/c< 



(21) 



where nrcc is the mean cumulative number of recombinations per 
hydrogen atom. 



6 CRITICAL POINTS IN THE 21 CM HISTORY 

Now that we have assembled the relevant physics, we are ready 
to construct some plausible histories. But first it is useful to iden- 
tify three "critical points" in the 21 cm history: the redshift Zh 
at which the IGM is heated above Tj, the redshift Zc at which 
the Wouthuysen-Field mechanism couples first Tg and Tk (when 
Xa = 1), and the redshift of reionization Zr. Their relative timing 
determines the structure of the 21 cm spectrum. 

The most important question is whether Zc occurs before the 
other two transitions: if not, the IGM will remain invisible. We can 
compute the collapse fraction (or more precisely A/coii) required 
to achieve Xc, = 1 from equation <18> . Substituting this into equa- 



tion <11> . integrated over time, we have the net X-ray heat input 

ATc at Zc-. 






0.08 fx 



fx,h /coll 

"02 
20 



A/coll 
3 



9690 1 0.6 0.023 

Nc Sc /roc fib/i^ 



(22) 



Note that ATc is independent of /,t because both the coupling and 
heating rates are proportional to the star formation rate (by assump- 
tion, of course). Interestingly, for fiducial (Pop II) parameters Zc 
should occur somewhat before Zh (see also the simple models of 
Chen & Miralda-Escude 2004; Flirata 2005). For Pop III stars, Na 
is a factor of two smaller and fx may be significantly larger than 
unity, in which case absorption would be much less important. Thus 
the presence or absence of a strong and extended absorption epoch 
offers a useful probe of high-redshift star formation. 

A similar estimate of the ionization fraction at Zc yields 



Xi,c 



0.05 



/esc iVion /coll 1 0.6 0.023 

1 + firoc Na A/coll Sa free ^bh? 

/_20_\2 
\l + z) ■ 



(23) 



For Pop II stars, Nion/Na ~ 0.45; thus even in the worst case 
(,fosc ~ 1 and nrcc = 0) coupling would become efficient dur- 
ing the initial stages of reionization. However, Pop III stars have 
much harder spectra and Nion/Na ~ 7-^' In principle, it is there- 
fore possible for Pop III stars to reionize the universe before Zc- 
This requires a small clumping factor, relatively late reionization 
(2 < 10), and /esc ~ 1- There is some theoretical motivation for 
such a large escape fraction in the earliest generation of galaxies 
iWhalen et all2004h . but maintaining that in larger galaxies seems 
unlikely (as might widespread and dominant Pop III star formation 
at 2 < 10). More conventional histories, in which /esc ^0.2 or 
2r ^ 12 fo r Pop III stars, predict s trong coupling throughout reion- 
ization (cf. lCiardi & Madaij2003l) . 

Finally, it is interesting to check whether we will see the IGM 
in absorption or emission during reionization. Combining equa- 
tions <11> and <2U . we find that the total heat input as a function of 



Xi IS 



AT{xi) / X. \ f fx,h 4800 10 

^ V 0.025 J r^ /esc iVion 1 



T., 



+ z 



(l + n,ec).(24) 



Thus, provided /x ^ 1 and Pop II stars dominate, it is safe to 
assume that the IGM is much warmer than the CMB throughout 
reionization. However, because Pop III stars have much higher ion- 
izing efficiencies, it is possible for X-ray heating to become signif- 
icant only when Xi > 0.2; in this case we would be able to see the 
early stages in absorption so long as fx remains small. This would 
be exciting in that it would add extra structure to the reionization 
signal, but it would also make isolating Xi{z) more difficult. 



^ Note that Fig. 1 of ICiardi et alj 120031) seems to suggest much smaller 
ratios for Pop III stars; this is because they report the number of photons per 
frequency interx'al. The frequency range between Lya and the Lyman limit 
is much smaller than that available to ionizing photons, which accounts for 
the difference. 
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Figure 1. Global IGM histories for Pop II stars. The solid curves take our fiducial parameters. The dot-dashed curves take /x = 0.2. The short- and long- 
dashed curves include photoheating feedback (see text), {a): CMB, gas kinetic, and spin temperatures (dotted, thin, and thick curves, respectively), {b): Ionized 
fraction, {c): Differential 21 cm brightness temperature against the CMB. In this panel, the two dotted lines show <5Ti, if shock heating is ignored. 



7 GLOBAL HISTORIES 

We will now construct some example histories using the prescrip- 
tions given in ^JSHS] Of course, given the uncertainties in each of 
the parameters, making firm predictions is all but impossible. But 
we can highlight the possible range of features. We will necessar- 
ily neglect a number of potentially important effects (such as metal 
enrichment and the evolution of the H2 cooling efficiency), but we 
hope that our histories are nevertheless representative of the possi- 
bilities. 



7.1 Pop II Stars 

We begin with a fiducial set of Pop II parameters. We take nimin = 
m4, A = 0.1, /esc = 0.1, jx = 1, A^ion = 4000, and iVc = 
9690. (Thus C = 40 for this model.) The solid curves in FigureQ 
show the resulting global history. In panel (a), the dotted curve is 
T-y, the thin curves are Tk, and the thick curves are Tg (note that 
they always interpolate between T^ and Tk). 

Figure^ shows the ionization history. Without feedback, it 
increases smoothly and rapidly over a redshift interval of Az ~ 
5. For these parameters, reionization is essentially complete at 
Zr ~ 7. Other values of ( move Zr back and forth but do not 
strongly affect its width in redshift space. Changing /* affects all 
three processes (nearly) equally: it simply moves the entire curve 
forward or backward in redshift. Note that the evolution slows 
somewhat at Xi > 0.8 and never quite reaches complete reioniza- 
tion; this is because the clumping factor C increases rapidly as 
Xi ^ 1 I Miralda-Escude et al. 2000). The dashed curves in this 
panel use photoheating feedback to slow reionization. Each as- 
sumes that the minimum virial temperature increases to Th = 
2 X 10^ K in i onized regions, nea r the upper limit of theoreti- 
cal expectations ( Diikstra et al. 2004). All other parameters are as- 
sumed to remain constant in photoheated regions. The two curves 



choose different prescriptions for the effective volume in which the 
feedback occurs; they bracket the possible effects of photoheating 
iFurlanetto & Loeb 2005).^ Clearly, Xi evolves monotonically even 
with this relatively powerful negative feedback. It does, however, 
slow the pace of reionization when Xi ^0.5, spreading it over a 
significantly longer time interval.* 

FigureQ: shows the 21 cm brightness temperature increment 
STt relative to the CMB. Here for convenience we have also la- 
beled the corresponding observed frequency i/ — !^2i/(l + z), 
where 1/21 ~ 1420 MHz. The signal clearly has interesting struc- 
ture. The most robust feature is the ~ 30 mK drop in STi, during 
reionization. For our fiducial parameters (and without feedback), 
that spans a frequency interval of ~ 35 MHz. In this model, re- 
combinations are slow (delaying reionization by only Az ~ 1) and 
Xi is essentially proportional to /coii- Thus, unless the ionizing effi- 
ciency increases throughout reionization because of some positive 
feedback mechanism (or if more massive galaxies drive reioniza- 
tion; Furlanetto et alj|2006l) . the signal is unlikely to have a much 
larger gradient during reionization. But with feedback, it can be 
significantly smaller: in that case, STt falls to zero over Az ~ 8 
(~ 60 MHz here). It is this gradient that might allow us to isolate 
the cosmological signal, so feedback will make its detection rather 
more difficult. 

A stronger feature appears near Zc and z^. At z> 20, the IGM 
is nearly invisible even though Tk <^ T^- Once the first galax- 
ies form, the Wouthuysen-Field effect becomes important and Ts 



^ In detail, the short-dashed curve assumes that the photoheated vol- 
ume is identical to the ionized volume, while the long-dashed curve cal- 
culates the photoheated volume without including recombinations. See 
IFurlanetto & Loebll200^ . 

* In fact, with these parameters feedback pushes Zr "^ 5, which is ruled out 
observationally. But that is of no consequence to us because we only wish 
to examine its qualitative effects on the 21 cm signal. 



© 0000 RAS, MNRAS 000, 000-000 



8 S. R. Furlanetto 



[It] 



103 



10^ 



10 

1 

0.8 
0.6 
0.4 
0.2 




- ' \' ' ' ^ 






















V, 


(a) 


'i; 


- 


; \ 


\ 




T, 


- 


^V" 






T, 


" 


— ^Os^ 


^N ^ 






— 














Xv ^ ^ 


.... ...— 




- 




^■^^>-'^ 


^^ 




- 




x^-^'^'C^''^ 






- 


- 


x^^ -. 


^ ^__^__, 


_^— " " 


- 


' 1 1 1 1 1 








: 


_ 1 1 1 1 1 








- 


: — \ 


~A 






; 


\ 


\ 




(b) 


- 


-_ \ 


\ 
\ 
\ 






: 


; \ 






I 


: \ 


\ 






- 


\ 


\ 






- 


- 


\^ \ 






- 




^\^.__ "^ .^^ 






- 


- , , , , 1 


1 1 1 1 1 1 


1 1 1 1 


1 1 1 


- 



50 



V (MHz) 
200 150 100 75 



10 



15 
z 



20 



25 



-50 



-100 



-150 



Ml 






- 


- 


.' 


(c) 


- 


y 








- 







- 


~ 


\\ / 




" 


- 


\^_j/ 




- 


1 1 1 


1 1 1 1 1 1 1 1 1 1 




- 



10 



15 



20 



25 



Figure 2. Global IGM histories for very massive Pop III stars. Panels are the same as in Fig.Q The soUd curve takes our fiducial Pop III parameters (without 
feedback). The long-dashed lines take /esc = 1. the short-dashed lines take fx = 5, and the dot-dashed fine (shown only in c) assumes /esc = 1 and 



drops rapidly. Because Na is large for Pop II stars, efficient cou- 
pling precedes the heating transition, allowing a relatively strong 
absorption signal {5Tb ~ —80 mK) to appear by z ~ 20. Ab- 
sorption persists until 2 ~ 15 (or ;/ ~ 90 MHz) before fading 
into emission. Thus, for these fiducial parameters, the universe is 
cold when coupling becomes strong, but it still heats up well before 
reionization begins in earnest. The absorption epoch can therefore 
teach us about the appearance and spectr al properties of the first 
sources of light in the universe JTozzi et a l. 2000; Sethi 2005), es- 
pecially the relative importance of Wouthuysen-Field coupling and 
X-ray heating. This era also produces much sharper features than 
reionization itself: the signal changes by -^ 110 mK over Az ~ 5 
(or /S.U ~ 20 MHz at such high redshifts). Thus, even though the 
foregrounds are also much larger at the lower frequencies, the early 
evolution may not be much harder to identify (especially if feed- 
back substantially slows reionization). 

Note that our model predicts a significantly stronger absorp- 
tion epoch than lChen & Miralda-Escude. 12004.1 : the difference lies 
in the assumed fx and in their assumption of a linearly increasing 
emissivity throughout reionization. In reality, structure formation 
proceeds exponentially at high redshifts, so the X-ray (and ioniz- 
ing) luminosities most likely decline much more rapidly with red- 
shift than assumed in their model. This permits a longer epoch of 
absorption. 

The other curves in Figure Q illustrate how the evolution 
changes with some of our input assumptio ns. First, the dot-dashe d 
curve assumes fx = 0.2 (near the value of lGlover & Branci2003h . 
This has no effect on Xi, but it delays the heating transition by 
Az ~ 3. As a result, Ts remains relatively close to T-, during 
the beginning of reionization, and ST^ never reaches its peak value 
in emission. On the other hand, the absorption trough is deeper, 
though it is still confined to z > 15. Also, we emphasize that the 
IGM remains visible throughout reionization, even when Tg < T^ . 

The two dotted curves in panel (c) show 5Tb if we ignore 



shock heating. In this case, all of the IGM material remains cold. 
At the high redshifts (z> 15) that are relevant here, shock heat- 
ing in filaments affects < 10% of the IGM, so its effects on 5Tb 
are modest. If weak shocks do occur outside of collapsing struc- 
tures, the absorption epoch will be much weaker. Moreover, if heat- 
ing is delayed until lower redshifts where structure formation has 
progressed much further, shocks should have a much larger effect 
(Gnedin & Shaver 2004; Furlanetto & Loeb2004J) . 

7.2 Population III Stars 

Figure |2| shows similar histories for Pop III stars. Here we take 
fiducial parameters (shown by the solid curves) m-mi-n = rn^, 
/* = 0.01, /esc = 0.1, fx = 1, and TVion = 30,000, yielding 
C, — 30; reionization ends at z ~ 7 in this case, ignoring feedback. 
The thermal history is similar to the Pop II case, except that the 
coupling and heating transitions are much more closely spaced in 
time because Pop III stars have relatively few soft-UV photons per 
high energy photon. As a result, Ts — > Tk occurs just about when 
Tk ~ T-y. Thus, if Pop III stars dominate, the absorption epoch will 
likely be considerably weaker - with a peak signal \5Tb\ < 50 mK. 
In general, with Pop III stars, we find that the time interval over 
which each quantity evolves is similar to the Pop II case, but the 
peak signals are up to a factor of two weaker because their spec- 
tra push the three transition points closer together. Thus detecting 
the global signal will be somewhat harder if massive Pop III stars 
dominate throughout reionization. 

The other curves use different input parameters. The long- 
dashed curves take /esc ~ 1 (which may be appropriate for the 
first stars, Whalen et al. 2004, although it is likely an overestimate 
in larger galaxies). In this case the universe is obviously ionized 
much earlier, but the thermal history remains the same. As a re- 
sult, reionization occurs while the IGM is still cold (see ^. This 
would complicate the interpretation of 5Tb{z) (because tempera- 
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Figure 3. Global IGM histories for a two-population model. Panels are the same as in Fig.0 Both sets of curves assume Pop III stars form in cold regions and 
Pop II stars form in hot regions (each with our fiducial parameters, except that /* = 0.1 for Pop III stars). 



ture fluctuations could be significant) but actually makes the reion- 
ization signal slightly stronger. The more important point is that the 
Wouthuysen-Field effect is already strong, so fluctuations will ap- 
pear throughout reionization: even though Nion/Na ~ 7 for these 
stars, coupling precedes reionization because each Lya photon is 
scattered many times. As suggested by equation <23t . aiTanging it 
so that Zc < Zr requires extremely unusual parameters, so it is al- 
most impossible to imagine that the IGM would be invisible during 
reionization. 

The short-dashed curves assume fx = 5. This strongly in- 
creases the heating rate, further decreasing the strength and dura- 
tion of the absorption epoch. On the other hand, this choice cleanly 
separates the reionization and absorption epochs and yields a brief 
period where Xi is small b ut the I GM s trongly emits. It gives a 
similar total heating rate to lHirata i2005h . who also predict weak 
absorption. 

Finally, the dot-dashed curve (shown only in panel c) takes 
/esc ~ 1 and fx = 5. This is, in many ways, a "worst case sce- 
nario" for studying reionization, because it compresses all of our 
critical points into a short time interval. Nevertheless, there remains 
interesting structure during and even before reionization. The dif- 
ference is that separating the effects of Ts, Tk, and Xi will be more 
difficult. 



reionization ipurlanetto & LoeblllOOSh - but it provides a simple 
way to implement exceptionally strong feedback. As in Figure Q 
the solid and dashed curves bracket the effectiveness of photoheat- 
ing feedback. In the fast case, there is a short period in which Xi 
overshoots because of the dramatic decline in ionizing efficiency 
between the Pop III and Pop II stars. The resulting recombination 
epoch is nevertheless shallow, and producing true "double reioniza- 
tion" scenarios requires pushing parameters in even more uncom- 
fortable directions. 

Nevertheless, Figure |3| shows that strong feedback can dra- 
matically slow the later stages of reionization. For obvious reasons, 
photoheating feedback becomes important when Xi ~ 0.5. Thus 
the gradient remains more or less unchanged during the first half of 
reionization (and the absorption epoch is completely unaffected), 
but it decreases strongly in the middle stages. Any potential "re- 
combination epoch" is also difficult to see, because Xi changes by 
only ~ 15% over Aa ~ 4 (yielding d&Tb/dv ~ 0.1 mK MHz"^). 
On the other hand, if feedback is fast, the late stages of reioniza- 
tion are completed entirely by the Pop II objects and so have as 
large a gradient as the single-population models. Thus strong feed- 
back could perhaps be inferred by observing two widely separated 
periods with large gradients. 



7.3 A Dual Population History 

In Figure |3| we exaggerate the strength of the photoheating feed- 
back transition by assuming that Pop III stars form in "cold" re- 
gions while Pop II stars form in hot regions. Both populations take 
our fiducial parameters except that we set /* =0.1 for Pop III stars. 
Thus the ionizing efficiency declines from C}^^ = 300 to (^^ — 40 
during photoheating, and TTimin simultaneously increases. There is 
no real physical justification for conflating the metallicity transi- 
tion with photoheating - in fact it is much more gradual than pho- 
toheating and requires fine-tinuing to synchronize enrichment and 



8 DISCUSSION 

In this paper, we have described the basic mechanisms relevant to 
computing the globally-averaged 21 cm signal from high redshifts. 
The basic structure of the background depends on three crucial 
transition points: when Ts -^ Tk (due to the Wouthuysen-Field 
effect), when Tk first exceeds T-y (most likely due to X-ray heat- 
ing, possibly together with weak shocks), and reionization. The rel- 
ative timing of these three transitions determines how much we can 
learn from the 21 cm transition. 

We have examined the range of possible signals and empha- 
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sized their qualitative features. If Pop II stars dominate tlie back- 
ground, and if tlieir X-ray emission (per unit star formation rate) 
is comparable to that in the local Universe [fx ~ 1), Lya cou- 
pling will become efficient before X-ray heating. This leads to a 
relatively strong absorption epoch (with a peak \5Tb\ ~ 100 mK) 
until X-rays cause Tk S> T!y. Once that occurs, STfy saturates 
until reionization begins - which is well after the heating tran- 
sition so long as /x ^ 1- Such histories contain long periods in 
which density fluctuations dominate the 21 cm signal, which is 
useful for studying cosmological parameters teow man et all2005t 
[McOuinn et al. 2005). This points to an important benefit of our 
calculation: it provides a basis for evaluating the expected ampli- 
tude of the 21 cm fluctuations during various phases of structure 
form ation. Cosmological measurements, spin temperature fl uctua- 
tions tearkana & Loeb'2005b';'Pritchard & Furlanetto'2005^, mini- 
halos (Iliev et al. 2002; Furlanetto & Oh 2006), and reionization it- 
self must all be placed in their proper contexts, and this is a first 
step in that direction. 

If very massive Pop III stars dominate reionization, we ex- 
pect a rather different history.^ Because these stars are so hot, the 
number of Lya photons per ionizing photon is much smaller. Thus 
the Lya coupling transition occurs much nearer to reionization, al- 
though the reionization epoch will still be visible in the 21 cm tran- 
sition. With these stellar populations, disentangling the effects of 
Ts and Xi may be more difficult. It would be even more compli- 
cated if the X-ray heating transition also occurs during reionization 
(which it could if we calibrate to the local X-ray luminosity-star 
formation rate relationship). On the other hand, they could easily 
have much larger X-ray luminosities, perhaps eliminating the ab- 
sorption epoch altogether 

A substantial quasar population (seeded by the remnants of 
massive Pop III stars, for example) could also eliminate the ab- 
sorption epoch, because such objects have large X-ray luminosi- 
ties. We have conservatively ignored this possibility because the 
kno wn po pulation of bright quasars contributes little to reioniza- 
tion JFan et aL2004) and does n ot affect the thermal history signif- 
icantly iVenkatesan et al.l|200ll), but this i s not to say that "mini- 
quasars" will be unimportant rOh 2001; Glover & Brand 20038; 
[Madau et al. 2004; Kuhlen & Madau 2005) - although they are un- 
likely to domin ate the ionizing photon budget liDiikstra et all2004t 
ISalvaterra et all2005.V These sources would only affect our predic- 
tions before X-ray heating from stars becomes significant - thus 
observing the absorption epoch (or lack thereof) could help to con- 
strain their significance. At lower redshifts, and especially dur- 
ing reionization, miniquasars are unimportant for the 21 cm back- 
ground because stars easily suffice to heat the IGM above T-y . 

Another uncertainty is our simple model for shock heating 
JFurlanetto & Loel3l l2004Y While it appears to provide a reason- 
able description of shocks produced as part of the cosmic web 
JKuhlen et alJl200a) . some simulations may imply that the low - 
density IGM could also host weak shocks JGnedin & ShaveI^2004^ . 
These are not included in the model and deserve further study. As 
with quasars, these would not affect our predictions during reion- 
ization, but they could shorten or even eliminate the absorption 
epoch. 

Clearly, measuring the mean 5Tb{z) would offer insight into 
the earliest generation of luminous sources {M adau et al. 1997; 
IShaver et alll999l;lTozzi et alJ2OO0l:lGnedin & Shavej2004 . How 



^ If metal-free stars have a more or less normal Salpeter IMF. they produce 
histories resembling those of our Pop II models. 
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Figure 4. Brightness temperature as a function of Xi in several different 
scenarios. We have removed the redshift dependence of the underlying sig- 
nal for clarity. The solid and short-dashed curves show our fiducial Pop II 
and Pop III cases, respectively. The dotted curve uses sets fx = 0.2 for 
Pop II stars, and the dot-dashed curve uses /esc = 1 and fx = 0.5 for 
Pop III stars. Finally, the long-dashed curve con'esponds to the solid curve 
in Fig.|3] 



usefu l migh t these measurements be for quantitative constraints? 
ISethil '200^ computed the sky-averaged 21 cm signal, including 
Lya coupling. X-ray heating, and ionization in a restricted range 
of models (including non-uniform X-ray heating but ignoring Lyn 
photons. Pop III stars, and feedback, for example). The qualita- 
tive features are similar to our Pop II models, although we predict 
a much wider range of possibilities given the uncertainties in the 
models. 

Treatin g the efficiency of each process as a free parameter, 
ISethil y003) considered how well these parameters could be mea- 
sured from STb{z). Because (with Pop II stars) heating generally 
occurs well before reionization, he found that Xi{z) could be well- 
constrained. We find considerably larger degeneracies. Figure |4| 
shows the brightness temperature as a function of Xi, with the ex- 
pected redshift dependence of STi, in a fully neutral universe re- 
moved. We show several different scenarios with a wide range of 
ionizing and X-ray heating efficiencies. To determine Xi precisely, 
we would hope that all these curves would overlap. Indeed they do 
approach each other for a;; > 0.5, but earlier in reionization the his- 
tories differ substantially. In most cases, the (1 — Ts/T^) factor 
has not yet saturated in this regime, either because of delayed Lya 
coupling (usually the case for Pop III stars) or relatively late X- 
ray heating (particularly in the dotted curve). Thus we predict that, 
in the absence of other constraints, Xi cannot be determined with 
certainty from the 2 1 cm background alone unless it is large. 

The most straightforward observation would be to measure 
5Tt{z) directly. Because this is an all-sky signal, single-dish mea- 
surements (even with a modest-sized telescope) can easily reach the 
required mK sensitivity. However, the muc h stronger synchrotron 
foregrounds (primarily from the Galaxy; see lFurlanetto et an200a 
for a review) nevertheless make such observations extremely diffi- 
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Figure 5. Gradient in the cosmological 21 cm signal for several models: 
our fiducial Pop II history (solid curve), our fiducial Pop III history (dashed 
curve), and a history with strong feedback (dot-dashed curve, corresponding 
to the solid curve in Fig.0. The small hitches in each curve occur at Tg ~ 



cult. As a rule of thumb, a "quiet" portion of the sky has a bright- 
ness temperature 



T, 



sky 



180 



180 MHz 



K. 



(25) 



Clearly, the amplitude of the cosmological signal will be impossi- 
ble to extract from the foregrounds. But the spectral structure of 
the 21 cm background imprints features on the (otherwise smooth) 
spectrum that may be separable. Figure|5|shows the expected gradi- 
ents in a few of our models. During reionization, we typically have 
d5Ti,/dz < 1 mK MHz^^, over three orders of magnitude smaller 
than the foregrounds. Without feedback, the gradients tend to in- 
crease throughout reionization; with feedback, they peak earlier on. 
The gradients can be much larger during the absorption era, partly 
because the signal can vary over a much wider range and partly 
because du/dz also decreases. Of course, the foreground gradients 
also increase significantly at these low frequencies (by about an or- 
der of magnitude at :^ « 90 MHz according to ea. l25> . so it is not 
clear which regime will actually be easier to constrai n. Separatin; 
the c osmological signal promises to be a challenge ( Sh aver et al 
ll999£.Gnedin & Shaver..2004.') and will be easier in histories with 
strong features (in which feedback does not slow reionization, for 
example). 

The li miting factors in dir ect measurements are likely to be 
systematics ( S haver et alJll999l) . Foregrounds vary across the sky, 
which couples to the frequency dependent beam and sidelobes 
(though the varia tion should be smooth and hence removable; 
lOh & Ma ck 2003 ). One strategy to avoid confusion with features 
in the foreground is to use multiple (wide) fields, in which fore- 
grounds are presumably uncorrelated but across which the cosmo- 
logical signal would be constant. Hydrogen recombination lines 
leave spectral features throughout the relevant frequency range, 
but they can be removed through follow-up observations with high 
spectral resolution. Finally, the biggest challenge will be calibrat- 



ing the telescope response to better than one part in 10'' over a wide 
frequency range. Internal loads of sufficient quality may be difficult 
to find. 

Fortunately, even if systematics do prove to be insurmount- 
able, there are other ways to extract the mean sienal. lCooravl i2005l) 
pointed out that scattering by a massive, low-redshift galaxy clus- 
ter can shift the spectrum in frequency space, so that differen- 
tial measurements between the cluster and background may re- 
veal the global spectrum. This avoids many systematics, but given 
the relatively small angular size of clusters may not measure the 
truly "global" signal. A second possibility is to use the 21 cm 
fluctuation amplitude to measure the mean backgrou nd, provided 
three-d imensional anis otropies can be measur ed t Barkana & Loebl 
l2005al) . For example, iMcOuinn et alj J2005h showed that, when 
density (and velocity) fluctuations dominate the signal, the overall 
signal amplitude can be extracted by combining 21 cm data with 
CMB data. (They phrased their constraint as Xi, but in more gen- 
eral models it corresponds to the mean brightness temperature.) 

Despite these difficulties, the global 21 cm background would 
be an extremely useful tool in understanding the overall evolution 
of the first sources of light. Together with the fluctuating signal 
(see lFurlanetto et all200d and references therein), the global back- 
ground could reveal many of the important properties of the "twi- 
light zone" of structure formation. 

I thank S. P. Oh for insightful comments on the manuscript, A. 
Loeb and J. Pritchard for helpful discussions, and R. Barkana for 
sharing the stellar spectra used in !J4| I also thank the Tapir group 
at Caltech for hospitality while some of this work was completed. 
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